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SUMMARY 


The Prandtl-Busemann small-perturbation method is utilized to 
obtain the flow of a compressible fluid, past an infinitely long wave- 
shaped wall. When the essential assumption for transonic flow (that 
all Mach numbers in the region of flow are, nearly unity) is introduced, 
the expression for the velocity potential takes 'the form of a power 
series in the transonic similarity parameter. On the basis of this 
form of the solution, an attempt is made to solve the nonlinear differ-' 
ential equation for transonic flow past the wavy wall. The analysis 
utilized exhibits clearly the difficulties inherent in nonlinear-flow 
problems. The investigation has nevertheless been rigorously carried 
to the point where the questioir of the existence or nonexistence of a 
mixed potential flow , free of disco^ntinuities can be settled by the 
behavior of a single power 'series Un the transonic similarity parameter. 
The calculation of the coefficient's of this dominant power series has 
been reduced to a routine computing problem by means of recvirsion for- 
mulas resiilting from the solution of the differential equation and the 
boundary condition at the wall. One of the interesting results of the 
analysis is the rigorous statement that the transonic similarity param- 
eter must be less than fo\ir-thirds. . - • 


INTRODUCTION 


, The present paper considers in detail a notion first expressed in 
reference 1 ; namely, that the Von KarmSn transonic similarity rule is 
implicitly contained in the potential-flow solution for high- subsonic 
flow past a prescribed body. The calculations are very much simplified 
by the’ choice of rectangular Cartesian coordinates as independent 
variables. For this reason the solid boundary chosen is a two- 
dimensional wavy wall of small amplitude extending to infinity in both 
the downstream and upstream directions. The problem is first treated 
by means of the Prandtl-Busemann iteration method for high-subsonic 
undisturbed speeds. On the basis of this solution an attempt is made 
to solve the nonlinear differential equation for transonic flow 
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(corresponding to the linear Tricomi equation in the hodograph plane) 
past a wavy wall. The analysis utilized exhibits clearly the diffi- 
culties inherent in nonlinear— flow problems. The investigation has 
nevertheless been rigorously carried to the point where the question 
of the existence or nonexistence of a mixed potential flow free of dis- 
continuities can be settled by the behavior of a single power series in 
the transonic similarity parameter. The calculation of the coefficients 
of this dominant power series has been reduced to a routine computing 
problem by means of recursion formulas resulting from the solution of 
the differential equation and the boundary condition at the wavy wall. 

The author wishes to acknowledge the invaluable aid and advice of 
Dr. A. Busemann of the Langley Laboratory during the writing of this 
paper, especially with regard to the final section "General Analysis." 

CALCULATION OF HIGH- SUBSONIC FLOW PAST A WAVY WALL 
BY MEANS OF PRANDTL- BUSEMANN ITERATION METHOD 


The fundamental nonlinear differential equation for the potential 
flow of a compressible fluid can be written as 


where 


(1) 


sL . 1 . 2_|J: M„^[l - (u^ * v^)] 

^00 

(2) 

and 



0 

X,Y 

velocity potential of flow / 

rectangular Cartesian coordinates in plane of flow 


u,v 

U 

c 

fluid velocity components along X- and Y-axis, respectively 
undisturbed stream velocity . 
local speed of sound 
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Coo speed of sound in undlstiurbed fluid 

Moo Mach number of undistvirbed stream 

7 ratio of specific heats at constant pressure and constant . 

volvune ; ' 

The quantities 0^ X, Y, u, and v are nondimensional with a char- 
acteristic length I as unit of length and the undisturbed stream 
velocity U as unit of velocity. The subscripts X and Y in equa- 
tion (1) denote partisQ. differentiation with respect to the designated 
variables. 

In order to obtain the Prandtl-Busemann iteration equations based 
on small perturbations of the undisturbed stream^ the assumption is made 
that the velocity potential 0 can be expanded in the form 


. 0 = X + 01 + 02 + 02 . .(3) 

For the purpose of defining and controlling the iteration procedure, the 
fiinction 0n+l and its derivatives' are regarded as small compared with 
the preceding approximation 0jj and its deriva,tives. Prom eqtiation (3) 
and the fact that for irrotatiohal or potential flow u = 0jf and 
V = 


and 


u = 1 + 


^IX + i^2X i^3X 


V - 01Y + 02Y + 03Y ■*■••• 

When these expressions for u arid v axe introduced- into equation (l), 
together with the corresponding expression for c^/c^^ given by equa- 
tion (2), and the powers and products of 0j^ and their derivatives are 
grouped according -to the assumptions of the small-perturbation method, 
the following iteration equatioris for the first three approximations 
result: 


^Ixx ^lyy “ ® 


(4) 
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^3xx * ^3yy 


+ 02yy = * ®)l»lxli*lxx 

(5) 

: ' 2M ^ X 

(1 +' cr) (03jj^^02x + ^Ix^^x) + 

w 

- : ■ 

[ 23^(1 + cj) - ^ 01x^ly^lxy 


(1 + ( 

j)|2p^(l + a) - ^0ixx^Lx^ 

i^^lixi^^ly^ + . 

i^lxy^2y ^ly5^2xy'\ 



where' x and y are new independent variables defined by the 
.transformation 


X = X 

y = 3Y 


(7) 


and 


= 1 - 


a = 


7+1 




Equation (4) is a Laplace equation and equations (5) and (6) are Poisson 
equations where the right-hand sides contain only previously determined 
quantities. These equations have been treated recently. in reference 2 
where the particular integrals of equations (5) and (6) are given in 
real form. These particular integrals have been utilized in obtaining 
the flow. over the wavy wall. Thus, the equation of the infinitely long 
wave-shaped wall (fig. l) is assumed to be - ‘ 


Y = a cos aX 


(8) 
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where ; . 

a amplitude of wave 



X . wave length . 

The reference element of length is conveniently chosen to be i 

The equation of the wavy wa.ll in nondlmensional form and in terms of 
the variables x and y then becomes 

y = Pe- cos X (9) 

where, if the thickness coefficient of the wavy wall is defined as 
t = then 6 = jtt. ' 

The expression for the velocity potential 0 , obtained with the 
aid of the particular integrals given in reference 2 , including terms 
of the order and satisfying the bo\indary conditions at infinity and 

at the solid wall to the same order, is as follows; 



0 = X + f e“^sin x + J 2 cM ^y + I2 + cr + (2 - ^e"^sin 2 x - 


l/£ 


8Ap 


-5 




+ 2a - llo^ - a(22 + 13 a) + 




■|a(6 + 5a) + (^ + 10a + 3a^)p^ e" 3 y + 4 a (2 + a)M^„Ve“^^sin x + 



6 M^^a ^3 + 2a + 


y sin 3x + 


( 10 ) 
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It is of interest to examine equation (10.) when the assumptions 
of transonic-flow theory are introduced. These assumptions are essen- 
tially that the imdistiarhed flow velocity differs only slightly from 
the speed of sound, that the velocity component normal to the oncoming 
flow is small compared with the speed of sound, and that the velocity 
component in the direction of the oncoming flow is. of the order of the 
critical velocity c*. If the undistiarhed stream is in the direction 
of the positive x-axis, then the velocity potential 0, referred to the 
critical velocity c*, can be written as (see reference l) 


0 = * + 7^(1 - M«2)f(x,y) (11) 


The second term on the right-hand side of this equation is the 
disturbance-velocity potential, and implies that terms involving powers 
of 1 - higher than the first are to be neglected. . The differ- 

ential equation satisfied by the function f(x,y) is obtained from-., 
equation (l) and takes the following simplified nonlinear form: 


^ ^ 
8y^ 8x 


( 12 ) 


The boundary conditions to be fulfilled by f(x,y) are as follows; 


Sx 




> 


'■ (at : y = «>) 


> 


— = -k sin X (at y = 0, ^ x ^ n) 

8y 


(13) 


where k is the transonic similarity parameter 


(7 + l)e 
(1 - 


m 
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Thus, it is seen that f is a function of x, y, and the transonic 
similarity parameter k only. In the limiting case, — > 1 and 

6 >0, k retains its meaning as a transonic similarity parameter. 

Note that in the Prandtl-Busemann- iteration procedure the order 
to which the boundary condition is applied at the surface of the profile 
is the same as the order to which the iteration has been taken. In the 
transonic case, however, the surface bovindary condition is always 
applied at the axis y = 0. in particular, this statement of the boimd- 
ary condition for transonic flow is a rigorous one in the limit as both 
the thickness coefficient and 1 - simultaneously approach zero but 

with the same distribution of slope that belongs to the family of pro- 
files being treated. A rather complete discussion of this point is 
contained in reference 3. 


Before application of the foregoing considerations to equation (lO), 
it should be noted that 0 in that equation is referred to the imdis- 
turbed stream velocity U. In order to Introduce c* as the reference 
velocity, both 0 and the right-hand side of equation (lO) must be. 
multiplied by u/c*. Now, from equation (2), with c = c* and 

u^ + v^ = - — , it follows that 
U2 




7+1 

2 





(15) 


where 


= JL 

c* 


If terms of only the first power of 1 - are retained, it follows 

from equation (l5) that . 


M* = 1 - 




(16) 


Then, multiplying the right-hand side of equation (lO) by this expres- 


sion for 


M*, replacing e/p by — equation (l4). 



. ^ o ■ - 

and neglecting all terms containing- 1 - ..to a power higher ;-than ,. . 

the first yields the following expression for 0: m • >' ■ 


0 = 


X + 


-^(l - "I"* ^ k3j,e7ysin x + 


(l + 2y)k^e"^sin 2x + (5 + %)sin x + 

+ 12y + 9y^)sin 3^k3e“3y^ + ... 


(17) 


Comparison of this equation with equation (ll) shows that the result 
obtained by means of the Prandtl-Busemann .iteratio.n method contains. . 
implicitly the form of th® solution required by transonic-flow theory. 
Moreover., the expression for the function f(x,y), namely, - ' 


f(x,y) = 


-X + 





,2y)k2e 


-2y; 


sin 2x + 


(5 + %)sin X + + 12y + 9y^)sin 3^k3e"3y + . . . 

(18) 

satisfies the nonlinear differential equation (12) and the boundary 
conditions, equation (13)., to the order k3. This fact suggests a 
solution for the flow over the wavy wall, in the neighborhood of Mach 
number Unity, obtained directly from equation (12). . , . . , 

. ■ SOLUTION OF EQUATION FOR TRANSONIC FLOW PAST A WAVY WALL 


Equation (l8) suggests the following' form for a solution of the 
nonlinear differential equation (12) for transonic flow subject to the 
boundary conditions stated in equation (13): . 
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f(x,y) = -X + ^ fjj sin nx (19) 

n=l 


where the fh’® functions of y only. This form for f(x,y) is 

substituted into equation (12) and repeated use made of the following 
identity: 


00 00 00 n 



When the coefficients of the separate harmonic terms sin nx are placed 
equal to zero,; the following system of nonlinear ordinary differential 
equations for fj^ results: • 


fn" - 

i n ^ (m + l)(m + n + Df^+ifm+n+l 
.■ m=0 

, (n = 1, 2, ... 00) ( 20 ) 

Note- that equation (20); with the right-hand side placed equal to 
zero, has a solution of the form e"^y. An iteration procediire is then 
set up in such a manner that the highest power of k and of e'^ 
appearing in the- functions f^ is equal to the number designating the 
iteration step. Thus,, the following setsof iteration -equations result: 

1st step, : 

•fl"'-fl=0 (21) 


= - X ZZ °- 

^ -m=0 


“)Vn-m - 




,.( 22 ) 


2d step. 
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3d step, 


4th step, 


fl" - 

f3" . 9f3 = .3f^f2 



(23) 


fg” - kf^ = - I fi2 - 3fif3 
- l6f4 = -6f;L^3 - kf2^ 


(24) 


5th step, 

fl" - fl . -fifj - 3fsf3 

/ 

fj" - 9fj . -3flf2 . 6fifi, 

- 25f^ = - 15f2^3 

'6th- step, . • . ■ • 


(25) 


■ V - ‘fs “ - 1 - 3 V3 - ® ■ 

f - I6f4 = -6f3_f3 - 4f22 - 

fg" - 36f6 = -I5fif5 - T ^3^ ' ^^^2^^^ 

and so forth. 



( 26 ) 


The right-hand sides of these equations are known functions con- 
structed from previously determined quantities in accord with the 
iteration procedure adopted. 
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These equations eire of the second order, non- homogeneous type with- 
const ant coefficients and are. readily integrated. The resulting ' : 
expressions for fj^, f2, fp f^, f^,, and . fg, with the houndary; 

conditions ( 13 ) talcen into account, are as follows: 


. ( 1 * il. >2 * 1861 aV.-y . r 

V 256 '^ * 3 X 2562 y 


? . t 8765 ;^2 -, 

256 9 X 256^ 


f + 33 k^'\y k3e”3y _ . 

\64 64 X 256 \ 576 X 256 ■ 72 X 256 


65 . : 7 • 

+ y + 


9 ' y 2 -^ 1 v3\k5e"5y + 

32 X 256 ^ 8 X 256 


(27) 


f, = 1^ . . ^^9 ^ ^ _ 2342 ^ ^4 , ( 1 , 

2 hr 72x256 432 X 2562 \Q 1024 


k'>) 3 k 2 e- 2 y . , 


24 X 256^ 


36 X 256 108 X 256^ 



815 ~ 

v2\;,r 

+ JL^ k2\v^ 

[256 

288 X 512 “ ■ \128 

8192 ^ rj 

i . 12245 

621? _ , 

- T9 „2' 

V 144 

X 2562 

’72 X .2562' . 

3 X ‘2562 : ' 

47 

v3 + 

1 

>e"^y +'• 

ro 

X 

256 

12 X 256 ^ J 



(28) 
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f3 = 


[ l ; 

[72 


. 23603 ^. J.p 

27 X 256^ 


+ 


\24 


72 X 256 J 


% * 


kV3y _ 


( 


1765 +• _151_.y + ll'i' 


3 X 256^ , 32 X 256 32 x 256 


‘y^ + ^ * * 


(29) 


^4 = 


Lli 

, . 390547 j.? 

+ /j^-+ 

32947 j 

[1536 

720 X '256^ ■ 

\384 

36 X 256^ 

/_3_^ 

617- j,2V^ 

■ + (— + 

11 ,J 

I128 

288 X 256 y 

V96 

24 X 256 


+ 


1 


( 


35251 , 2465 455 

^ Jr 


90 X 256^ 576 X 256 . 96 X 256 




119 , y3 + . ... ... >V6e-6y ^ 

i»8 X 256 ^ -W X 256 ^ / 


(30) 


•p. .= /-I— + -I- y + y2 + y3 + ^ yAk^e'^y + 

5 I38IK) 768 ^.. ■1536 384 24 X 256 . 


(31) 




91 91 

+ — ^ 7'y 




^6-: \J^^54 X 2048 72 x 256 ^ 20^48’^ 1024 ^^^ 


y4 + , 9 , y5V6g-6y •+ 
2048 ^ 5120. ^y . 


(32) 
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IS 


Note that the functions f^, f 2 ) aJid. f^ include the terms of equa- 
tion ( 18 ), obtained from equation (lO) by allowing the Mach number to 
approach unity. ^ ^ ! 

Equations (27) to (32) may be considered to be essentially the 
nonlinear solution for the flow past a wavy wall of small amplitude for 
stream Mach numbers in the neighborhood of unity in the form of a power 
series in the transonic, similarity parameter k. Moreover, this, solu- 
tion is identical with the one obtained by means of the- Prandtl-Busemann 
iteration equations when Mach number unity is approached. 

CALCULATION OF LOCAL MACH NUMBER, CRITICAL SIMILARITY 
PARAMETER) AND PRESSURE COEFFICIENT 
General Formulas 

From equation (ll), when all terms containing 1 - Moo^ to a power 

higher than the first are neglected, the expression for the fluid 
velocity referred to the critical speed of sound c* is 



1 + 


— (1 

+ 1\ 




( 33 ) 


Now, from equation (2), the relation between . q and the local Mach 
number M is 



7+1 

;2 




1 + 



• ( 34 ) 


In the transonic approximation, the difference of any Mach nvimber in the 
field of flow from \inity is considered a small quantity. If terms of 
only the first power .of 1 - M^ are retained, equation' ,(34) yields 



1 + 


2 

7+1 



+ . . • 


( 35 ) 
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Hence , from equation . ( 33 ) ^ 


1 _ m2 = -(1 - m„2)|£ ■ . (36) 

This equation provides a means of calculating the critical value of the 
transonic similarity parameter, that is, the value of the parameter .k 
for which M = 1 at. the point of maximum velocity on the boundary. 


Calculation of the Critical Value of k 

For the family of wavy walls of small amplitude (including the 
limiting case of vanishingly small amplitude) with M = 1 at x = 0 , 
and y = 0, equation (36) yields the following relation for the deter- 
mination of the critical value of k: 



y=0 


(37) 


■' By means of equation (19) together with the expressions for f-]^ 
to fg ' given by equations (27) to (32), equation (37) yields the 

following power series, exact to seven terms, for the determination of 
the critical value of k for the family of wavy profiles:' 




^k3 

384 


337 

9216 




' '4o43 
256.x, '576 


k5 


359381 

, r 

270 X 256" 


*,= 1 

(38) 


The procedure adopted in order to estimate the critical value of k is 
as follows: ' Froifl'equktion ( 38 )' the- value Of k' cto be found for 2, 3, 
5] 6, ‘and 7 terms. 'These values' of ‘ k ; are, ' respectively, 1, O.899O, 
0.8644, ■ 0;‘85047 0.8424, and 0'.8377i The last two valiaes indicate 'the ' 
approach to the asymptotic value of k, that is, the value of k when 
the mimber of terms in equation (38) is infinite. If the values of k 
approach smoothly to the asymptote,, the estimated critical value 
k = 0.8377 is very nearly correct. ; 7 " 
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Suppose nov that "both sides of equation (36) are divided "by 
Q7 + 1)0 Then, since 

k = 

(i-M.2)3/2 


equation ( 36 ) can be written as 



v-2/3 


( 39 ) 


The right-hand side of , equation ( 39 ) is a function of x, y, and the 
parameter k only and is. characteristic of the entire family of boundary 
profiles. For the family of wavy walls with the critical value assigned 
to k, equation ( 39 ) evaluated at the wall becomes 


[(7 + l)0^/3 


1.1253(1 0.8536 cos X - 0.0989 cos 2 x - 0.0299 cos 3 x - 

0.0122 cos ifx - 0.0038 cos 5x - O.OOl? cos 6x) (4o) 


Table I lists the values of — ; — - — j- for values of x between 

■[(7 + 1)0^'^ 

-jt and jt and figure 2 shows the corresponding curve. 


Calculation of the -Pressure, Coefficient ; 

Bernoulli' s theorem for a compressible fluid assimies the . following 
differential form along a streamline:; 



. l 6 
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where 

p pressure in fluid 

p density of fluid ' . ' ' 

q speed of fluid 

and the subscript 00 denotes the quantity in the \mdisturbed fluid. 

Now, from equation (2) the relation between the nondimensional speed q/U 
and the. local Mach number M is given by 



1 .+ 


1 + 





i^2) 


In accordance with the assimption of transonic -flow theory that all Mach 
numbers in the flow differ only slightly from unity, equation (42) 
becomes , 



2 

7+1 





(43.) 


where powers of 1 - Mq, and 1 - M higher than the first have been 
neglected. Then with the pressure coefficient defined as 




P - Pcp i 

I pj^' 


equation (4l) takes the -following form' for transonic flow: 


dC. 






After integration. 





7 + iL 



(44) 
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IT 


or, with the aid of equation (36), 




hk - 8 


(^ 5 ) 


The right-hand side of this equation can he considered to be the first 

. - • ■ , . ■ . • 2 • ■ 

term in a power- series development of Cp^j^ in 1 - , In particular 

when the local Mach n\jmber- first attains ' unity, then , ^ = 0 at ' x = 0, 


y = 0 and equation (45) becomes 


8x 




cr 


Tl(^ ■ “4 


,(46) 


a result valid in the transonic range only. Again, if the thickness 
coefficient approaches zero as Moo,— — >1 and - M = 1 at x = 0, y = 0, 


.dC- 


then equation (46) shows. that .the slope ' 




d(l - uj) 


: at the critical 


value Moo = M, 


cr . 


1 is a constant - 


-, independent of the par- 


7+1 

ticular family of profiles treated. ,This result is valid whether the 
approach to Mach'huinber unity is made from the subsonic or supersonic 
region.' — • ' ^ . 

• (7 + 1 -)^'^^ ’ • 

If both sides Of equation (45) are m\iltiplied by : , then 

■ , ■■ . 2/3 : . 


(7 + 1 ) 
2/3 


1/3 


S,Moo = 


-2/3 


t • i) 


(47) 


where the right-hand side depends only on x, y,.. and -the parameter k, 
and is characteristic of the ’family of bo-undary profiles treated. For 
the family of wavy walls with the critical- value of- k-' chosen, equa- 
tion' (4 t) takes the following form at the solid' siarface (y = 0) : 
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(7 + 




= -2.2507(0.8536 cos X + 0.0989 cos 2 x + 0.0299 cos 3 x + 
0.0122 cos 4 x + 0.0038 cos 5x + 0.0017 cos 6x) ( 48 ) 


(7 + 1)^/^ , 

Table I lists the values of — - — « for values of x between 

g2/3 ' 

-rt and n and figure 2 shows the corresponding curve. ' 


GENERAL ANALYSIS 


An examination of the expressions "for f]L 
tlons (27) to (32) shows that the general form 


to fg given by equa- 
of f^ is 


p=0 


-(2p+n)y 


n-1 

2p+n_2 00 


q=l 


y : a” ? w 


n+2r 


r=p 


(h = 1 , 2 y . . . 00) (49) 

where> if p = 0, the upper limit of q is n-1, and, if p ^ 0, the 

upper limit of q is 2 p + n - 2 . The four-labeled coefficients P 
' q r 

are real numbers calculated from recursion formulas obtained from the 
system of differential equations (20) and the boundary condition at the 
surface of the wavy wall. The boundary Condition at y = 00 is auto- 
matically satisfied by the form of f^; whereas the boundary condition 


at the wall takes the form 


(^n')y=o 0 , (n 1) ^ 

. Inserting the expression for f^^ given by equation ( 49 ) into equa- 
tion (50) yields immediately the following results: 


1 0 
0 0 


= 1 
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and 


(2P + h)Aq 1= fZ 

■ .P= 0 , 1 . 


(n — Ij 2 f . « « 00) (51) 


where, if n = 1, the lower limit of p on the right-hand side is unity, 
and, if n / 1 , the lower limit of p is zero. Also, if n = 1 , the 
upper limit r of p goes from 1 to 00, and, if n / 1,. r goes - 
from 0 to 00. 


In order to find the recxarsion formulas for the coefficients A^ 

q r' 

the expression f or . f ^ given by equation ( 49 ) .'is substituted into the 
system of differential equations (. 20 ). The calculation is facilitated 
by the introduction of the following notations; 


n-1 

2 p+n — 2 00 _ 

Aj. r ( 52 ) 

^ r=p 

and ' ' ' 


A^<q, ^ X / ,(53) 

where the quantities A™^’^ arise from the multiplication of the two 

00 00 

infinite series ^ a“ e"^Py and >_ A^ e'^P^. Note' that the quan- . ' 

p=o ^ ■ p=o ; 

titles are symmetric with respect' to' the upper labels m and n. 

Then • .. . ' 
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When this expression for is, suhstituted into equation. (20.) ,. and- 

repeated use is made of the identity 


n 


Bm,n - EE ®ni,n 

m=0 n=in n=0 m=0 


the exponential terms in y can be eliminated, and the following' reciir- 
Sion formulas result: 



and 

l4-p(p + n)A^ - 2(2p + n)^A^y + (aJ)" 

p 

1 ^ — / ■ \i2m,m.n+m 

- n 2_ A p-m 

m=l 

(P = 1, 2, 


(p = 0; n = 2, 3, . . . ~) (55) 



h-1 

E 

m=l 


m(n 


myA’^^ 


n^m~ ■ • . •' -■ 

P 


. 00; n = 1, 2, , . . 00) ( 56 ) 


Note that these recursion formulas still contain powers of y and k. 
Both y and k . can be eliminated for given values of p and recursion 

formulas containing only the coefficients Ag ^ thus obtained. For 

example, consider equation (55) for which p = 0. By repeated use of 
the identity 



t t 

3=0 a=3 • 
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the following recursion formula is obtained: ■ 


-2n(q + l)Aq+i ° + (q + l)(q '+ 2)6^"^ A^2 ? = 



n-2-t 


r ( 

i=n _ -h 


m + 


l)(n 


m 


- 1) 


s=0 




0 


0 

r-s 


(■^ ” 3 } • • • Q.' — O^ 'X^' 2^ • « • (n — 2)j 

r = 0, 1, 2, ... 00) (57) 


where 


(q = n - 2) 
(q / n - 2) 


n-2 0 

’ q ■ 1 


A number of interesting relations can be easily obtained from this recur- 
sion formiila. Thus, for q = n - 2 and r = 0, 


2(n - 1)A 
n 


n 

-1 


0 

0 


1 

4 


m=0 


- m 


1) (m + 

n-2-m 


0 ^m+1 0 

0 m 0 


Then 



1 

8 


0 ^ 2^ 

0 ” 48 X 128 




32 


.6 0 

A5 0 = 


5120 


7 (58) ■ 


4 0 
3 0 


96 


From these numerical values the general form is found to be 


A_“ ° 

'^i-l 0 


n' 


n-2 


niV 


n-1 
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Similarly, for q = n - 2 and r .= 1 , 


, n 0 _ n^-2 ^1 0 . 


In a corresponding manner, more complicated expressions can be obtained 
„ . n 0 , ■ „ „ • 

for Aj^_ 2_ ^ when r = 2, 3,, ■. ... 

Note that in the expression for the disturbance-potential function 

00 

f(x,y), infinite series of the type ^ y^"^e"’^yk’^'*’^stn nx A^^^ ® 

occur. The ratio of the (n + l)st and nth terms is 


kye 


-y sin(n + l)x 


sin nx 


(^n+1 0 
^ n '1 

k\^' 

n-1 1 


In order for: the series. to converge, the limit of this ratio as 'n — > «> 
must be less than unity. Thus, if the maximum value of ye"^^= is ' 
inserted and x = 0 , Cauchy's ratio test yields 


k e T ,1 

— -i-<l or k <4 
e 4 


( 59 ) 


Other infinite series occur in the expression for f(x,y) which 
diverge for values of the transonic similarity parameter k considerably 
less than 4 . Thus, consider the recursion equation ( 57 ) and take q = 0 
and r = 0 . Then 


- 2 nA 


;n 0 


1 0 


+ 26 


n-2 

0 


A 


n. 0 1 

2 0 = ■ T ^ 


n-2 

y~ (m + l)(n 

m=0 


m - l)A 


m+1 0 
0 0 


,n-m-l 

A 0 


0 

0 
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and 


a2 0 ^ 1 

10-8 


6a| q = 2a| 0 + 3A§ 0 

§ + 3a3 0 + 2(a§ 0)' 


4a^ ° = a^ ° 
1 0 2 


lOA: 


5 0 


^5 0 


if 0 




i“i 0 = - ¥(^ g)‘ 


• • « • 


Also from equation (57) > 

1 .3 0 2 0 

4 A 2 0 = Al O' 

8 A 2 0 3A^ § + Uf g A§ g + 3 a| § 

20 a| g ^ 6 a 5 g . i5Af g a 3 g . 15 a 3 g 4 g . ioaJ g 

8 a| g = 2 A 3 g + 8 Ai g A^ g + 9 a 3 g Aq g + 8 aJ g Ag g + 5 A^ g 


and 


30a5 ° = 13aJ ° ^ 15a| 0 A= ° ♦ ioaJ ° + i5a| ° a3 ° 

12a| ^ 4a6 0 , 9 ^ 0 ,3 0 , 8 A,^ 0 ,2 0 , 5^5 0 , 8,2 0 ,4 0 , |(,3 ^ 

48A^ ° = 20A« ^ 24A^ ° a| ^ 15a’ ° ^ 27a| ° 4 I * ^>>4 S 4 ? 
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The following relations are obtained from the supplementary equations (51) 


2A^ ^ = A^ ^ 
0 0 10 


5 0 _ 5 0 
'0 0 " 1 0 


3a3 0 = a3 0 
0 10 


0 a6 0 

S 0 =•*10 


ta** ° . A** ° 
■^0 0 1,0 


From these relations, together with theipnes listed in equation ( 58 ), 


A? ° = i 
^ 0 8 


a5 0 =,-L 

1 0 768 


A? ° = ^ 

^ 0 2k 


pf> 0 ^ 7 X 13 
10 72 X 256 


A.^ 0 = i 

*1 0 . 381 , 


and lead to the general rule 


K 


n 0 
1 0 


{3° - ^ 


.(n = 2, 3, 


(60) 


From equation ( 5 l), nA^ q = A^ q, and, -therefore, , 


' n o {3n - 5I " 

^ 0 = ^ ,n-/ (n = 2 , 3 , . . . 
. nnl4^ ^ 


(61) 
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. wbgre, by, definition, 



.4 X 7 X 10 X 13 X 


3n - 5 


In the expression for the local Mach nximber in the field of flow given 
by equation (36), the expression 6f/6x occ\irs. This expression, with 
X = 0, contains infinite series of the type 


F 


00 


nk“e‘ 




0 

0 


and 


G 


7" nk^e-^^A? 


0 

0 


( 62 ) 


(63) 


If the maximum value of e~^y(= l) and the expressions given by equa- 
tions (60) and (61) are inserted, the Cauchy ratio test shows that 

k 

both series converge for k < — . In particular, the series expression 

for F evaluated at the surface (y = o) can be expressed in closed 
form. Thus, 



The graph of 
point at k = 


F 

4 

3 


against k is a semicub ical parabola with the cusp 
and F = 2. With the restriction that the transonic 


similarity 'parameter k be positive and that one and only one value 
of k correspond to a given value of F, the permissible values of k 
and F are confined to the part of the parabola lying between the 
origin (0,0) and the apex (4/3,2). The power-series expression for G 
evaluated at the surface can also be expressed in closed form; namely, 
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/ X 1 n-1 

(n - 1)I4 


- i 


- k 


(65) 


This expression, together with the one for F, shows clearly that the 
parameter k cannot he equal to but must be less than four- thirds, 

A close examination of the recursion formulas (55) and (56) dis- 
closes the important fact that each one of the manifold of power series 
in k that appear in the functions fj^ can he expressed in terms of 
the members of a single dominant set of power series. This dominant 
set consists of one power series in k from each f^^ namely^ the one 

multiplied by e only. According to equation (^9), the members of 
the dominant set of power series are given by 


Sn = A? r (n = 1, 2, . . . co) (66) 

r=0 

Several examples are now given to illustrate this important obser- 
vation. Consider the series that belong to the set 



n+2r 


(n = 2, 3, 


oo) 


From the recursion formulas. 


oo 


2 0 2r+2 _ 1 o 2 

^1 r ^ " 8 


C *1 r ^ * I 


r=0 


w 

r 

r=0 


k 0 2v+k _JLs.^ + 1 q2c .3„_ ^1.2 

1 ^ ^ ^ 16 + 4 S1S3 + - S2 
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Consider now the set 


00 

E 

r=0 


•n 0 n+2r 
^2 r ^ 


(n = 3, 


Then 


> a 3 0 k2r+3 = i_ 

• 2 r 32 


r=0 


00 

5 


4 0 , 2r+4 _ 1 q 4 . 1 o 2q 

2 ^ ^ ' I28 4 ^2 


Consider the set 


00 

E 

r=l 


,n 1 , n+2r 
Aq r ^ 


(n = i; 2 ^ . • . 00) 


From the recizrsion formula (56), 


00 

1 ; 


,1 1 . 2rtl _ 3,3 1 , , 


.2 1 , 2r+2 
Aq r 4 


17 q^_J_q 2 „ _lqq 

12 X 256 ^ 96 ^ ^ 413 




3 1 ,^2r+3 ^ 3^5 4^ S,^S^ - Sn^So - — StSo^ - ^ S-.S) 


256 ' 


2 ^ ^1 ^3 - 3 i ^ 1^2 -.8 
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Consider the set 


£ 

r=l 


A° ^ 
1 r 


(n — 2^ • • • 00) 


Then 


00 

r=l 


^1 1 j^2r+l ^ Q 3 


1 r 


64 


00 


2 1 ,.2r+2 3 g ^ _ 1 g 2g 


Ai ; k 


256 ^ 8 


1 °2 


^31 i, 2 r +3 _ 15 o. 5 57 o 3 cj 9 . 2 „ 3 o cj 2 

^ ^2 - ^ ^1 S3 - Ta Si®2 


r=l 


2048 ■ 512 1 2 - ^ -1 -3 - 16 "1^2 


other examples are the following: 




a1 2 k2r+l . 1 s 5 

3 r 2048 -‘- 


,1 2 2r+l 
*2 r ^ 


3 o 5 3 o 3 c 

8192 " 256 2 


^ a 1 2 k2r+l = _J_ S.5 1 _ s. 3 s 2 - -i SnSp2 - i/s.3s 

^- 7 ^ 1 r 8192 ^ 1536 ^ 2 T 1 2 64 2 


r=2 


16 


64 


1 2 ij2r+l ^ 61 




256 X 576 ^ 384 


1^ + ^ S£3s 2 - i Sj^S2^ '+- 


^ ^ SgS^ 
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The recognition of the existence of a dominant set of power series 
in k represents a major reduction in the complexity of the present 
problem. Thus, consider the array of infinite power series contained 
in equation (66): 


Si = 4 ° ° .3 . ° k? . ° ic-' . . . . 

S 5 . Ag 0 ^5 , ^ 0 J ^ ^ 0 ^9 , ^ 0 ^ _ 


( 67 ) 


An examination of equations (2?) to ( 32 ) shows that the coeffi- 
cients Aq ^ of the series appear to be positive and monotonically 

decreasing. The series formed from the first column on the right-hand 
side of equation ( 67 ) therefore dominates the series formed from sue- 

00 

ceeding columns. Moreover, the first-column series ^ ’ A? 0 has a ■ 

n=I 

radius of convergence k = (see equation ( 61 )). The radii of con-, 

vergence of the other columnar series therefore are either equal to or 
greater than four-thirds. Similarly, an examination of the series Sn 
shows that the coefficients in each column on the right-hand side of 
equation ( 67 .) also seem to form positive and monotonically decreasing - 
sequences. -This behavior means that is the dominant series of .the 

set Sjj and, in fact, of the aggregate of power series in k in the 
expression for the disturbance potential f(x,y). 

Consider now- the series consisting of the first- terms of the- odd- 
labeled series S^; S^, ' S^, . . ., that is. 
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n=0 


a2h+ 1 0 v2n+l 
^ 00 ^ 


( 68 ) 


According to the theory of power series (and it can be easily verified)^ 
the radius of convergence of this power series is still k = ^. 

Now, a comparison of corresponding terms with the dominant series 
shows that 


Ai ° = ^19 ^ 0 = Ji. 

^ ^ 72 X 256 ^ ° 72 


and 


1 0 
^0 2 


33^15 ^ .5 0 _ 7 

1*32 X 256 ^ ^ ° ■ 15 X 256 


Thus, if in general A^ ^ >A^q^ q, then the radius of convergence of 

the dominant series can be less than the radius of convergence of 

the comparison series given by expression (68) and therefore may con- 
ceivably be equal to the critical value . kj,^. = O.8377. Moreover, it- 

would then follow that the original Prandtl-Busemann small-pertxirbation 
method is valid for purely subsonic flows only. This conclusion would 
not invalidate other approaches to the transonic-flow problem 
(reference 4). 

Unfortunately, the coefficients A^ ® do not conform to any 

apparent or superficial law, but perhaps a carefta study of the recur- 
sion formulas (55) and (56) and the supplementary or boundary' rela- 
tions (51). will yield a rigorous proof of the foregoing statements. 
Otherwise, it remains to calculate a reasonable number of the coef- 
ficients Aq For this purpose the development of complete recin-sion 

formulas similar to' equation (57) for the required values of p is 
worth while. Thus, for p,= 1, the recursion formula is 
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Mn + 1 )a“ - 2(n + 2)(q + 1)6° A^J^ + 


n . n 1 


(q * l)(q + 2)6“-^''' . 


- ^ “^q (“ + 1)(“ - ” - 1) Z!^-t ° r-i+1 * 

t=0 m=q-t L s=0 


(--"■) Vt ^1 A”r r°s - 5 + 1) aJ ° 


^n+1 0 

q r-s 


(jn — 1^ 2^ • • 0 00 j ^”0^ 1^ 2^ • « • 

r = 0, 1, 2, ... 00) ( 69 ) 


where 


.n-l,n ^ 0 

q . 1 


(q = n - 1 or n) 
(q ^ n - 1 or n) 


Finally, it may be of interest to give the general formvila for the 
Mach number distribution at the surface of the wavy wall. Thus, by 
means of equations ( 36 ) and (^9), 



2m)cos(n 


m 

2m)x 

p=0 



.n-2m 


P 

m 


(70) 


where 


denotes the integral part of 


n/2. 
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TABLE I.- VALUES OF 


Qr t 


(7 + 1)^/^ 

AND — -L n at 

2/3 P'”cr 


THE SUBFACE FOR THE FAMILY OF WAVE- SHAPED WALLS 


I - 

\j7. ^ 

1.996929 

1.99^505 

1.987188 

1.957^30- 

1.906695 
1,83^90 
1.7i^3^33 
1.634666 
1.510815 
1.373602 
1. 300524 
1.224843 
1. 146916 
1.067146 
.903656 ■ 
.737110 
.569667 
.404761 . 
.250223 
.119935 . 

.031487 

.007977 -• 
. 000007 ■' ■ 
.007977 , 

.03148? 

.119935 
.250223 
.404761 
.569667 ■ 
.737110 
.903656 
1.067146 '■ 
1. 146916 
1.224343 
1.300524 
1.373602 
1.510815- 
1.634666 

1.743433 

1.834890 

1.906695 

1.957430 

1.987188 

1.996929 


(y !)• 


1.743190 ■ 
1.738342 

1. 723708 
1.664191 
1.562722 

1. 419111 

1.236197 • 

1.018664 
.770962 
. 496536 
-350379 
. 199018 
.043163 
-.116375 
-.443357 
-.776449 
-1.111335 
-1.44ii45 
-1.750223 
-2.010799 
-2.187694 
-2.234715 
-2.250654 
-2.234715 
-2.187694 
-2.010799 
-1.750223 
-1. 441145 
-1.111335 

-.776449 

-.443357 

-.116375 

. 043163 
. 199018 
.350379 
. 496536 
. 770962 
1.018664 

1.236197 

1. 419111 

1.562722 
1. 664191 

1.723708 
1.743190 








Abstract 


The simplified nonlinear differential equation for 
transonic flow past a wavy wall is solved by the method 
of integration in series. The solution has been carried 
to the point where the question of the existence or non- 
existence of a mixed potential flow can be answered by 
the behavior of a single power series in the transonic 
similarity parameter. The calculation of the coeffi- 
cients of this dominant power series has been reduced 
to a routine computing problem by means of recursion 
fomulas resulting from the solution of the differential 
equation and ^the boundary condition at the surface of the 
wavy wall. 
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